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Abstract 

In this paper we introduce new methods to prove the finite cyclicity of some graphics 
through a triple nilpotent point of saddle or elliptic type surrounding a center. After ap¬ 
plying a blow-up of the family, yielding a singular 3-dimensional foliation, this amounts to 
proving the finite cyclicity of a family of limit periodic sets of the foliation. The bound¬ 
ary limit periodic sets of these families were the most challenging, but the new methods 
are quite general for treating such graphics. We apply these techniques to prove the finite 
cyclicity of the graphic ( 1 ^ 4 ), which is part of the program started in 1994 by Dumortier, 
Roussarie and Rousseau (and called DRR program) to show that there exists a uniform up¬ 
per bound for the number of limit cycles of a planar quadratic vector field. We also prove 
fhe finife cyclicify of fhe boundary limif periodic sefs in all graphics buf one fhrough a 
friple nilpofenf poinf af infinify of saddle, ellipfic or degenerate fype (wifh a line of zeros) 
and surrounding a center, namely fhe graphics (/g^), and {Dl 2 b)- 


1 Introduction 

This paper is part of a long term program to prove the finiteness part of Hilbert’s 16th problem 
for quadratie veetor fields, sometimes written H(2) < 00 , namely the existenee of a uniform 
bound for the number of limit eyeles of quadratie veetor fields. The DRR program (see paper 
[2]) reduees this problem to proving that 121 graphies (limit periodie sets) have finite eyelieity 
inside quadratie veetor fields, and the long term program is to prove the finite eyelieity of all 
these graphics. 

This program has been an opportunity to develop new more sophisticated methods for ana¬ 
lyzing the finiteness of the number of limit cycles bifurcating from graphics in generic families 
of C°° vector fields, in analytic families of vector fields, and in finite-parameter families of 
polynomial vector fields. In this paper, we focus on some graphics in the latter case: graph¬ 
ics through a nilpotent point and surrounding a center inside quadratic systems. The general 

*This research was supported by NSERC in Canada. 
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Figure 1 : The graphics (J^^4) 




(C) (i/?3) (d) {Dhb) 
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method is to use the Bautin trick, namely transforming a proof of finite cyclicity of a generic 
graphic into a proof of finite cyclicity of a graphic surrounding a center. This is possible in 
quadratic systems since the center conditions are well known: indeed all graphics through a 
nilpotent point and surrounding a center occur in the stratum of reversible systems. The sys¬ 
tems of this stratum are symmetric with respect to an axis, and are also Darboux integrable with 
an invariant line and an invariant conic. In practice, the Bautin trick consists in dividing a dis¬ 
placement map 1/ in a center ideal, i.e. in writing it as a finite sum of “generalized monomials” 
times non vanishing functions of the form 

n 

V{z) = ^aimi{l + hi{z)), ( 1 . 1 ) 

i=l 

where each a* belongs to the center ideal in parameter space, rrii is a generalized monomial in 
2: and hi{z) = o(l) behaves well under derivation. 

To compute the displacement map, we write it as a difference of compositions of regular 
transitions and Dulac maps near the singular points. The Dulac maps are calculated in nor¬ 
malizing coordinates for a family unfolding the vector field. In this paper, we develop some 
general additional methods, which allow to prove the finite cyclicity of the graphic (144) (Fig¬ 
ure 1 (a)). In particular, for the unfolding of this graphic, it is very helpful to be able to claim 
that all regular transitions are the identity in the center case. This is possible if we exploit the 
fact that the centers occur when the system is symmetric, and if we choose cleverly the sections 
on which the different transition maps are defined. Also, in the center case, the Dulac maps 
have a simple form since the system is Darboux integrable. 

The methods can be summarized as follows. 

• We highlight that the change to normalizing coordinates in the neighborhood of the 
singular points on the blow-up locus can be done by an operator. This allows preserving 
the symmetry in the center case when changing to normalizing coordinates. 

• We introduce a uniform way of calculating the two types of Dulac maps when entering 
the blow-up through a much shorter proof than the one given in [8]. 

• Although each Dulac map is not C^, we can divide in the center ideal its difference to the 
corresponding Dulac map in the integrable case. 
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• The method of the blow-up of the family allows reducing the proof of finite cyclicity of 
the graphic to the proof that a certain number of limit periodic sets have finite cyclicity. 
These limit periodic sets are defined in the blown-up space. The ones obtained in blowing 
up a nilpotent saddle are shown in Table 2 . For all of them but one (the boundary limit 
periodic set), we can reduce the displacement map to a 1-dimensional map, the number of 
zeros of which can be bounded by the Bautin trick and a derivation-division algorithm on 
a map of type ( 1 . 1 ). The boundary limit periodic set is more challenging, since we need 
to work with a 2-dimensional displacement map, the zeros of which we must study along 
the leaves of an invariant foliation coming from the blow-up. We introduce a generalized 
derivation operator, which allows performing a derivation-division algorithm on functions 
of the type 

n 

V{r,p) = + hi{r,p)), ( 1 . 2 ) 

i=l 

where hi are C*^-functions on monomials and rrii are generalized monomials in r, p (see 
definitions in Appendix II). During this process, we have to take into account that rp = 
Cst. 

We have a partial result for every graphic, but one (namely (^^^4)), through a triple point at 
infinity: 

Theorem 1.1. Let us consider the graphics and {Dl 2 b) through a triple point 

at infinity (see Figure 1). Then for any of them, the boundary periodic limit set obtained in the 
blowing up has a finite cyclicity. 

Theorem 1.1 is not sufficient to prove that the given graphic has a finite cyclicity inside the 
family of quadratic vector fields. The reason is that, beside the boundary limit periodic set, 
other limit periodic sets (see for instance Table 2 for (144)) are obtained in the blowing up and, 
as explained above, we have to prove that each of them has also a finite cyclicity. We present 
here a complete result for the first graphic: 

Theorem 1.2. The graphic ( 144 ) has a finite cyclicity inside the family of quadratic vector fields. 

As for the finite cyclicity of the other graphics {Ilf), {Hlf and {Dl2b), we intend to address 
the problem in the next future. The finite cyclicity of (Hff) should be straightforward with ar¬ 
guments identical to those used for (144). It will be done simultaneously with the corresponding 
generic graphic {Hf2)- Some of the limit periodic sets to be studied for {Ilf will involve four 
Dulac maps of second type. For these limit periodic sets, it is not possible to reduce the study 
of the cyclicity to a single equation. Hence, new methods will need to be adapted to treat the 
center case, when the periodic solutions correspond to a system of two equations in the four 
variables ri, pi,r2, p2, with r4p4 = z/4 and r2P2 = 1^2- As for the graphic {Dl2b), some of the 
limit periodic sets to be studied involve four Dulac maps of second type, two of them through 
the semi-hyperbolic points P4 and P2 on the blown-up sphere. 

The techniques developed in this paper can be adapted for studying the boundary limit 
periodic sets of graphics of the DRR program through a nilpotent finite singular point. The 
only new difficulty in that case is to show that the three parameters of the leading terms in the 
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displacement map do indeed generate the center ideal. We also hope to adapt them to study 
the boundary graphic of the hemieyele there, the additional difficulty is the two semi- 

hyperbolie points along the equator. 

Proofs of Theorems 1.1 and 1.2 are given in Seetion 3 and Appendix II, where the detailed 
eomputations of cyelieity are found in Theorems 5.8, 5.12 and 5.13. Theorem 4.1 in Appendix 
I, gives a statement about normal form for 3-dimensional hyperbolie saddle points in a way 
adapted to this paper. Theorem 4.5 of the same appendix gives a new proof for Dulae transitions 
near these saddle points, shorter than the one given in [8]. Precise properties for the speeific 
unfoldings dedueed from the quadratie family are proved in Appendix III. These properties of 
some parameter funetions are needed to obtain the results of finite cyelieity. 


2 Preliminaries 

2.1 Normal form for the unfolding of a nilpotent triple point of saddle or 
elliptic type 

We eonsider graphies through one singular point, whieh is a triple nilpotent point of saddle or 
elliptic type. A germ of veetor field in the neighborhood of such a point has the form 

X = y 

( 2 . 1 ) 

y = ±x^ + bxy + yx'^y + yO{x^) + 0{y‘^). 

The saddle ease eorresponds to the plus sign, and the elliptie ease to the minus sign with 
|6| > 2a/2. In the elliptie ease, we limit ourselves here to the ease |6| > 2 a/2, whieh eor¬ 
responds geometrieally to a nilpotent point with hyperbolie points on the divisor of the quasi- 
homogeneous blow-up. 

The unfolding of such points has been studied by Dumortier, Roussarie and Sotomayor, [4], 
including a normal form for the unfolding of the family. A different normal form has been used 
in [8] for studying the finite eyelieity of generic graphics through such singular points, when 
we limit ourselves to |6| > 2^2 in the elliptie case. This normal form is particularly suitable 
for applieations in quadratic vector fields, where there is always an invariant line through a 
nilpotent point of multiplieity 3. 


A germ of C°° veetor field in the neighborhood of a nilpotent point of multiplieity 3 of 
saddle or elliptic type ean be brought by an analytie ehange of eoordinates to the form 

X = y + ax^, 

( 2 . 2 ) 

y = y{x -f- yx + o{x ) -f 0{y)). 

This requires an additional ehange of variable and sealing eompared to what has been done in 
[8]. The point is a nilpotent saddle when a < 0 and a nilpotent elliptie point when a > 0 (see 
Figure 2). The ease |6| = 2 a/ 2 eorresponds to a = 
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(a) Saddle case 


(b) Elliptic case 


Figure 2: The different topological types 


For a 7^ a generic unfolding depending on a multi-parameter A = (/xi, /i2, /xs, /i) has the 
form 


X = y + a{X)x^ + fi2, ^2 3 ) 

y = fii + fi3y + x^hi(a;, e) + y{x + t]x^ + x^h2{x, A)) + y‘^Q{x, y, A), 

where hi{x,X) = 0(|A|). Moreover, hi,h 2 ,Q are C°° functions, and Q can be chosen of 
arbitrarily high order in A. 

2.2 Finite cyclicity of a graphic 

Definition 2.1. A graphic F of a vector field Xq , i.e. a union of trajectories and singular points, 
has finite cyclicity inside a family Xx if there exists G £ > 0 and 5 > 0 such that any 
vector field Xx with | A| < 5 has at most N periodic solutions at a Hausdorjf distance less than 
e from T. If a graphic has a finite cyclicity, its cyclicity is the minimum of such numbers N. 

This means that when studying the finite cyclicity of a graphic F, we need to find a uniform 
bound for the number of periodic solutions that can appear from it, for all values of the multi¬ 
parameter in a small neighborhood W of the origin. Typically we need to find a uniform bound 
for the number of fixed points of the Poincare return map or, equivalently, for the number of 
zeros of some displacement map between two transversal sections to the graphic. With graphics 
containing a nilpotent singular point there is no way to make a uniform treatment for all A G W, 
and we typically cover PF by a finite number of sectors, on each of which we give a uniform 
bound. The method for doing this is the blow-up of the family, which was first introduced in 
[6], and next applied to slow-fast systems in [1]. 

2.3 Blow-up of the family 

We take the neighborhood of the origin in parameter-space of the form x [0, z/q) x 17, where 
t/ is a neighborhood of 0 in /x-space and we make the change of parameters 

(/Xi, /X2, /Xs) = (z/^/Xi, Z/^7X2, ZZ/Xg), (2.4) 
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where M = (/ii, /r 2 , /Xg) G and v G [0, uq). 

Note that is compact. Hence, to give an argument of finite cyclicity for the graphic F, 
it suffices to find a neighborhood of each M = ( 71 ^, 712 , 713 ) G inside a corresponding 
1^0 > 0 and a corresponding U on which we can give a bound for the number of limit cycles. In 
our study, we will consider special values Oq of a. It is important to note that a(A) depends on 
A, and hence that a — oq is in some sense a parameter in itself. 

The way to handle this program is to do a blow-up of the family. For this, we introduce the 
weighted blow-up of the singular point (0, 0, 0) of the three-dimensional family of vector fields, 
obtained by adding the equation z> = 0 to (2.3). The blow-up transformation is given by 

(x,?/, z/) = (rx,r^y,rp), (2.5) 

with r > 0 and {x,y,p) G After dividing by r the transformed vector field, we get a 
family of C°° vector fields Xa, depending on the parameters A = {a — oq, M, p). The foliation 
{u = rp = Cst} is invariant under the flow. The leaves {rp = u}, with u > 0, are regular 
two-dimensional manifolds, while the critical locus {rp = 0} is stratified and contains the two 
strata (see Figure 3): 

• X M+ is the blow-up of Xq (for A = 0); 

• Dji = [of -\-y‘^ + p^ = l\p>A], for any J1 G 

2.4 Limit periodic sets in the blown-up family 

The vector field Xa has singular points on r = p = 0. For a 7 ^ there will be four distinct 
singular points (occuring in two pairs) corresponding toy = 0 (for Pi and P 2 ) and y = 

(for P 3 and P 4 ): see Figure 3. Their eigenvalues appear in Table 2.4. 




Figure 3: The stratified set [rp = 0} in the blow-up 


We will study the finite cyclicity of a graphic V joining a pair of opposite points Pi and 
Pj+i in X, with i = 1 or i = 3. We consider a particular value Aq = (ao,Mo,po)- Here 
is the strategy for finding an upper bound for the number of limit cycles that appear for A in 
a neighborhood of Aq. We determine the phase portrait of the family rescaling (2.6) on D-p: 
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r 

P 

y 

Pi 

—a 

a 

-(l-2a) 

P2 

a 

—a 

(1 - 2a) 

P 3 

1/2 

-1/2 

-(l-2a) 

Pa 

-1/2 

1/2 

(1 - 2a) 


Table 1: The eigenvalues at Pj (i = 1,2,3,4) 


this allows determining limit periodic sets T, whieh are formed by the union of T with a finite 
number of trajeetories and singular points on joining Pi and Pj+i, so that their orientation 
be eompatible with that of T. The limit periodie sets to be studied appear in Table 2 for the 
saddle ease. They eome from studying the phase portrait of the family rescaling 


X = y + ax^ + p,2-i 

y = di+Pfy + xy, 


( 2 . 6 ) 


obtained by putting p = 1 and r = 0. It then suffiees to show that eaeh limit periodie set has 
finite eyelieity, i.e. to show the existenee of an upper bound for the number of periodie solutions 
of XA for A in a small neighborhood of Aq. 

2.5 Proving the finite cyclicity of a limit periodic set 

Typieally, the kind of argument we will use for proving the finite eyelieity of a limit periodie 
set is the following: we look for the zeroes of a displaeement map between two seetions. The 
seetions are 2-dimensional but, beeause of the invariant foliation, the problem ean be redueed to 
a 1-dimensional problem and the eonelusion follows by, either an iteration of Rolle’s theorem, 
or its generalization, namely a derivation-division argument. The teehnique ean be adapted 
to non generie graphies oeeurring inside integrable systems: the proof in the generie ease is 
transformed into a proof for the eorresponding graphie, using some adequate division of the 
eoeffieients of the displaeement map in the ideal of eonditions for integrability. 

To eompute the displaeement map, we deeompose the related transition maps between see¬ 
tions into eompositions of Dulae maps in the neighborhood of the singular points and regular 

transitions elsewhere. 


2.6 Dulac maps 

The Dulae maps are the transition maps in the neighborhood of a singular point on r = p = 0. 
They are eomputed when the system is in normal form. The normalizing theorem is The¬ 
orem 4.1 of Appendix I. There, it is proved that the normal form is obtained by a normaliz¬ 
ing operator J\f, a erueial property for this paper. The theorem establishes the existenee of a 
parameter-depending loeal ehange of eoordinates of elass bringing the blow-up of (2.3) in 

—N 

the neighborhood of one of the points P* into the normal form (up to f i—)■ —t) written in 
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Table 2: Convex limit periodie sets of hh-type for a graphie with a nilpotent saddle. 


normal form eoordinates (F, r, p) (provided that the eigenvalue in r has a sign opposite to the 
two other eigenvalues). Using Table 2.4, we take a = 2(1 — 2a) near cxo = 2(1 — 2ao) for P 3 
and P 4 when oq < |, and a = near uo = for Pi and P 2 when a > ^. The normal 

form X^ is given by 

1. If (To ^ Q : 


{ f = r, 

tr 

Y = -(a + (^A(z^))F. 


(2.7) 


8 






















2. If (Jo = I G Q, with {p, q) = I when q I : 

{ r = r, 

P_= -A _ _ (2.8) 

F = - (a + ipA{i^))Y + rPY'^)Y + ffVA{i^), 

with pa = 0 when (Jq ^ N ( g 7 ^ 1 ). 

The funetions if a, ^a,Va are polynomials of degree < K{k) increasing with k, with smooth 
coefficients in A and 0) = 0. 

We introduce the “compensator” function a), also denoted cUaiO^ defined by 

f aAO 

a;(e,a)=cuFO= (2.9) 

1 ^-lnF a = 0 . 

We propose in Appendix I a new computation of the Dulac maps previously studied in [ 8 ] . 
There are two types of Dulac transitions. The first type of transition map goes from a section 
{r = ro} to a section {p = po}> or the other way around. This type of transition typically 
behaves as an affine map, which is a very strong contraction or dilatation. The study of the 
number of zeroes of a displacement involving only Dulac maps of the first type is reduced to 
the study of the number of zeroes of a 1 -dimensional map. 

The second type of Dulac map is concerned with a transition map from a section {Y = Yq} 
to, either a section {r = tq}, or a section {p = Po}. We take z/q = roPo- 


2.6.1 First type of Dulac map 


Theorem 2.2. We consider the Dulac map from the section {p = po} to the section {r = tq}, 
both parametrized by (Y, u). Let 

d = d{a, u) = a + ipa{i^) 


and 

a = a{a, v) = d{a, 0) — (Jq. 

The Y-component of the transition map Da has the following expression: 

1. Ifao ^ Q : 

Da{Y,u)= (-YY. ( 2 . 10 ) 

2. Ifao = I G Q with (p, g) = 1 when (Jq ^ N : 

= PA(i^)Po(^) + (^) (y + (l>A{y,J^)), (2.11) 

with Pa as in (2.8). In particular, = 0 when ao ^ N. 

The function family (j)A in (2.11) is of order | In z/|) and for any inte¬ 
ger I > 2, is of class C^~‘^ in {Y, v^Puj a j , v, p, a). 
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2.6.2 Second type of Dulac map 


Theorem 2.3. Wfe consider the Dulac map from the section {Y = To}, parametrized by (r, p) 
to a section {r = tq} parameterized by {Y, a). It has the form (r, p) h- {DA{r, p), u), with its 
Y-component, {DA{r,p), given by: 

1. Ifao ^ Q : 

D^(r,p) = (-)Vo. ( 2 . 12 ) 

Vro/ 

2 . Ifao = I e Q with (p, g) = 1 when cto ^ N : 

+ (^) (Yo + (j)A{r, p)'^, (2.13) 

with Pa CIS in (2.8) (pa = 0 when ao ^ N). 

The function family (f)A in (2.13) is of order | In r|) and, for any inte¬ 
ger I > 2 , is of class C^~‘^ in a), p, p, a). 


3 Applications to quadratic systems 

3.1 Quadratic systems with a nilpotent singular point at infinity 

Theorem 3.1. A quadratic system with a triple singularity point of saddle or elliptic type at 
infinity and a finite singular point of center type can be brought to the form 

( X = -y + Box"^, 

\y =x-^xy, 

with Bq > 0. For Bq 1, the full 5-parameter unfolding inside quadratic systems is given with 
B = Bq + pq inside the family 

( X = -y + Bx^ + p 2 y‘^ + (pi + Bpfj x 
\ y =x + xy + psy'^ + (1 - 2B)p^y. 

For Bq = 1, the full 5-parameter unfolding inside quadratic systems is rather given with B = 
I + Po inside the family 


( X = -p + (1 + pq)x‘^ + p2y‘^ + P5X 

\y = X + {p4 + P5)x^ + xy + psy'^. 

The parameter p 2 (resp. ps) corresponds to a nonzero multiple of the parameter p 2 (resp. pY 
in the blow-up of the family at the singular point. There is no parameter pi in this family since 
the connection along the equator is fixed. 

Moreover for (3.1) we have: 
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1. Bq > 1 for a nilpotent saddle; 

i?o = I corresponds to a = —\ in (2.2) (6 = 0 in (2.1)). 

2. Bq < Ifor an elliptic point; the elliptic point is of larger codimension, type 1 (the singular 
points in the blow-up coallesce by pairs) if Bq = \ (corresponding to a = \ in (2.2), i.e., 
h = 2^/2 in (2.\)). 

3. The system (3.2) has an invariant line y = —1 if — (1 — 2B)p^ = 0. 

4. If p 2 = /^3 = /^4 = 0, the system (3.2) has an invariant parabola 


2B-1 


y = 


+ {2B - l)p^x - — + {2B - l)pl. 


(3.4) 


2 - ■ V -/r— 25 

The parabola y = \ is invariant for system (3.3) when po = p 2 = Ps = Pi = 0. 

5. The integrability condition is ps = Pi = P 5 = 0, for which we have the following 
graphics with return map 

• B>1: (II,). 

. 1 < B < 1/ (/i), 

. 0 < B < i; (Bfj), 

. B = O; (Hf,), 

. B = 1; (Dl2t). 

6. The value of “a ” in the corresponding normal form (2.3) is a = 1 — B, and the parameters 
P 2 and p 3 correspond to p 2 and p^ up to a nonzero constant. 

Proof. We can suppose that the nilpotent singular point at infinity is located on the y-axis, the 
other singular point at infinity on the x-axis, and the focus or center at the origin. Then the 
system can be brought to the form 


X = 5iqx + 5oiy + 52ox‘^ + 5iixy, 
y = 7ioa; + 7oi2/ +7113^2/+ 7022/^- 


(3.5) 


Localizing the system (3.5) at the singular point at infinity on y-axis by u = ^, w = ^, we 
have 


V = (^11 - 702)^^ - + (^20 - 7 n ) v ^ + (^10 - 7ol)^^^^.' - 


(3.6) 


For the singular point (0, 0) of system (3.6) to be nilpotent, we should have = 702 = 0. The 
point is triple if 711 f 0. 

We want the finite singular point to be a center, which corresponds in this case to the system 
being reversible with respect to a line. Because of our choice of singular points at infinity this 
line can only be the y-axis. Then 610 = 701 = 0. 

By a rescaling and still using the original coordinates (x, y), we obtain the system (3.1). 
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The change of coordinates W = —w + {Bq — l)r>^ brings the system (3.6) into the equivalent 
form 

j V =W 

[w = {Bo-l)V^ + {2Bo-3)VW+ o{V^)+ o{VW). 

The classification of the nilpotent singularity at infinity follows. 

A general unfolding preserving the singular point at the origin (which is simple) is of the 
form (after scaling of x, y, and t) 


f X = —y + Bx^ + rriiQX + muxy + mo 2 |/^ 
\ y = x + xy + noiy + n 2 ox^ + r^ 022 /^ 


with B close to Bq. We use a change of variable (X, Y) = {x + (ly, C 2 X + y) for small (^i, (, 2 - 
The terms in XY in the expression of X and the term in in the expression of Y vanish 
precisely when 

f (25 — l)Ci ~ ''^11(1 + C1C2) + 2(21x102 + 2Cin,o2(Ci + C2) — C1C2 = 0, 

1 (^ - 1)C2 + (1 + Ci)no 2 - drill + Clmo 2 = 0, 


which can be solved for (Ci,C2) by the implicit function theorem except for 5o = 1. When 
Bq = 1, we replace the second equation by the vanishing of the term in Y in in the expression 
of Y, namely 

Cl + C2 - noi + mioCiC2 = 0 . 

Again, we get a system that can be solved for (Ci, C2) by the implicit function theorem. □ 


3.2 Finite cyciicity of the boundary iimit periodic sets of {I 14 ), (Imj) and 

(Dh,) 

Notation 3.2. In the whole paper, * denotes a nonzero constant, which may depend on some 
parameters. 

Theorem 3.3. The boundary limit periodic sets of {lob) o-nd {Dl 2 b) (see Figures 1 (a), (b) 

and (d) and 4) have finite cyclicity. 


Proof. The finite cyclicity of the boundary limit periodic set is studied inside the family (3.2) 
when 5o 7 ^ 1, and we will discuss later the adjustment when 5o = 1. 

Choice of parameters. We take as parameters 

M = ifis, Pa, P5,fi2, Bo - 1) = {Mc,fi2,Bo - 1), (3.9) 

with {ji 2 iV“i) ^ Si and {Bo — ^,Pa, Pb) in a small ball. The parameters 

Me = {P3, Pa, Pb) (3.10) 

unfold the integrable situation. We let Ic be the ideal of germs of -functions of the parameters 

generated by {p 3 ,Pa,P 5 }- 
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Figure 4: The boundary graphic through P3 and P4 and the four sections Ej and Hi, i = 3,4, in 
the normalizing coordinates. 


Notation 3.4. 1. The symbol Op{Mc) refers to a function in the parameter M belonging to 

the ideal Ic. 

2. The symbol Og{Mc) refers to a function of(X, M) which belongs to the ideal generated 
by Ic inside the space of functions of (X, M). Depending on the limit periodic set, we 
could have X = x^,, where is the normalizing coordinate near P 3 , or X = (r, p). 

The displacement map. It is better to consider the chart y = 1 in the blow-up. We take 
normalizing charts in the neighborhood of P3 and P4. As discussed above, these normalizing 
charts can be chosen symmetric one to the other under the center conditions. The normalizing 
coordinates are (r, p, Tj) near P,. We consider sections Ej = {xi = Xq} and 11* = {r = ro} 
in the normalizing charts. The sections E* are parameterized by (r, p), and the sections 11* by 

{Xi, u). 

Let V = P 4 o S' — T o P 3 be the displacement map from E 3 to 114 : T and P 3 follow the 
flow forward, while S and P4 follow the flow backwards. 

Let us first give the proof when (T*(0) ^ Q. The Dulac maps are defined from sections 
E* = {xi = Xo} to sections 11* = {r = ro}, with Xq and tq fixed. Then the Dulac maps P* 
have the form 

P*(r,p) = (C'*(MK%rp). (3.11) 

We can choose Xq and tq so that C*(0) = 1, i.e. XgrQ= 1, and Cs{M) = C 4 {M) under the 
center conditions. 

The map T has the form 

T(T 3 , u) = (P(T 3 , u), u). (3.12) 
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Because of the symmetry of the sections, then H = id under the center conditions. 

The planes r = 0 and p = 0 are invariant under the map S, which hence has the form 

S{r,p) = {rF{r,p),pF-\r,p)), (3.13) 

with F of class C^, since u = rp is invariant. Moreover, it is known from [8] that F(0, 0) = 1 
when the sections Sj are symmetric. 

The displacement map then has the form 

A(r,p) = {C^{My^F^\r,p) - H ,v) . (3.14) 

Let V (r, p) be the first component of A. Then periodic solutions correspond to zeroes of V. 

We now need to compute F and F. 

Computation of H. 

The map H is in (Ts, v). It has the form 

ny, z/) = Ts + £o(M) + ei{M)x 3 + 0{xI)0g{Mc), (3.15) 

with eo{M) = Op{Mc),ei{M) = Op{Mc). 

For /r2 = = /^4 = 0, the system (3.2) has the invariant parabola (3.4). The term p 4 ^x in x 

is without contact, which yields that 

£o{M) = *p4{l + 0{M)) +0{p3) + 0{p5)0{M) = *p 4{1 + 0{M)) +0{py +0{p3)0{M), 

(3.16) 

where * denotes a nonzero constant. Lemma 6.1 in Appendix II shows that the same is true for 
(3.3). Let us again take P 2 = = 0. The divergence is then {2B + l)x + (1 — B)p 3 . 

Proposition 6.2 in the Appendix II shows that 

£'i(iff) = *^3^1 + 0[M)) + OiyP^v) + 0[p4). (3.17) 

The center ideal. The equations (3.16) and (3.17) imply that we can take {eo, as gener¬ 

ators of the center ideal Ic- 

Computation of F. The function F has the form: 

F{r, p) = l + w(l + 0{p)) + OyOciMc). (3.18) 

Indeed, it is proved in Lemma 6.3 in the Appendix that the second derivative of pF{0, p) is a 
nonzero multiple of /Ig. Moreover, the blown-up vector field is integrable on r = 0 for /Ig = 0. 

Writing the displacement as a finite sum of terms. We need grouping all terms of the dis¬ 
placement map into a finite sum of the form (1.2). We will see that three terms are sufficient 
and show that 

V{r,p) = -eo(M)(l -f Kir^p)) - C 3 {M)ei{M)F^{l + /ii(r,p)) -f *Jiy^p{l + h 2 {r,p)). 

(3.19) 

We now explain how to group the different terms. 
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Notation 3.5. The symbol O(r^) used in the sequel, is for an unspecified 6 > 0, which may vary 
from one formula to the other. 

Let us first consider the terms coming from H o D^. Remember that H is the identity when 
we have a center. Moreover, the map H really takes place in the initial {x, y)-plane, where the 
center ideal is generated by /^s}- Hence, the higher order terms of H o are of the 

form 

(eo{M)ko{r, p) + ei{M)ki{r, p) + P3k2ir, p)) ■ 

The first two terms contribute to ho{r, p) and hfir, p), as contributions of order 0{r^). As for 
the third term, we use the fact that = rpjl^. Hence it contributes to h 2 {r, p), also as a term of 
order 0{r^). The term will be later grouped with the corresponding term 

coming from D4 o S. 

Let us now consider the other terms coming from D 4 o S'(r, p) = CfiM)r'^'^F{r,pY'^. Again 
we use that F is the identity when there is a center, i.e. all its terms are divisible in the ideal 
Ic- One of them is the term *]I^r'^^p coming from the term *P3P of F. As mentioned above, all 
higher order terms r‘^'‘o(p) have coefficients divisible by JL^. Also, all terms in r'^'^pOfr) can be 
distributed in ho, hi and /i 2 , as terms of order 0(r‘^). Hence, we only need to consider the pure 
terms in It suffices to show that all such terms can be divided in This comes 

from the fact that the computation of the pure terms in r can be done in the plane p = 0, and 
that the system restricted to this plane does not contain any term in /Z3. Since 

d 4 -d 3 = iyOp{Mc)f{u) = rpOp{Mc)f{u), (3.20) 

with / of class C^, we can replace everywhere by d^, up to terms of order distributed 

in ho, hi and h 2 . 

We are left with the terms We write this as 

= (CoiM) - CfiM))r^^ + CfiM){r^^ - 

= {CoiM) - + CfiMfido - diy^u}{r,d 3 - 0 ^ 4 ). 

(3.21) 

The difference C 3 {M) — CfiM) is Xorg “^^(l — Using (3.20), the two terms can be 

decomposed in sums of terms contributing to ho, hi, /i2, as terms of order O(r^). 

Finite cyclicity in the case cxo irrational. The displacement map 1/ in (3.19) is a special case 
of a universal family 

ao(l + ho{r,p)) + air^^{l + hi{r,p)) + a 2 r'^^ p{l + h 2 {r,p)), (3.22) 

with ho, hi of order 0{r^) and /i 2 is of order 0(p) + 0(r‘^). Using that these three functions are 
of order o(l), we show in Theorem 5.8 below that this family has at most two small zeros along 
any curve rp = Cst for r, p < 5 for some small 5. This implies that, either V has at most two 
small zeros, or V is identically zero, in which case we have a center. 

Adjustment of the proof when ctq = | with g > 1. The adjustments are minimal. Indeed, the 
formula of the Dulac map is more complicated: 

A(r,p) = (r^(C',(M) + 0(r,p)),rp), (3.23) 
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with (j){r,p) as in Theorem 2.3. Hence p) produces in V new terms of order 0{r^), dis¬ 
tributed in ho,hi,h 2 . 

Adjustement of the proof when = p. Here the first component of Di{r, p) has an additional 
term of the form 

( T _ 

ro 

All higher order terms can be distributed in hQ,hi,h 2 and we need only consider the term 

E = o S — {1 + ei(M))K3 = (^4 o S' — K 4 ) + E with E = p) — (1 -1- ei{M))K-i{r, p). 



1. We consider first the term ^4 o S' — K4. Let (3 = — p. We have that 


K/i{rE) — K4{r) = 




\ro J 


Uji 


ro/J 


~G^ 


Let us consider G{r, p): 


Since 

V ro / 

we obtain that 

, . - 1 / r \ 

G{r,p) = - 


) — OJ13 i — 1 = 1 — 


^^0 


r \~3 E ^ — 1 


ro 


P Vro 


^ +{E^-l)uJ^) = 


ro 




E 0 - 1 


/S 


r \-0 


ro 


- - h 


ro 


i.e. G{r,p) = and then Ki{rE) — ^(r) = 

As F = 1 + */i3p(l + pg{p)) + rOciMc), we have that 

F/3 _ 1 


/3 


= *p3p{l + pp(p)) + rOciMc), 


and then that 


Ki{rE) - Ki{r) = -p4Z/Pr^(*p3p(l + pp(p)) + rOciMc)). 

The term rOciMc)) gives contributions of order 0{r^) in ho,hi,h 2 . Next, the term 
*P3p(l + pg{p)) gives the contribution — * piuP~^p{l + pg{p)) in /i2. If p > 2, this 
term is also of order 0(rp), and it is of order 0(p) if p = 1. 

2. We consider now: 

E = pP 


{P4{jy) -p3{jy){l + ei{M))r''^uj y-^as-p 
+ P4(i^) {r""^ a3 - p^ 

+PA{jyy^ ( a; ( —, ^3 - p ) - a; ( —, ^4 - p 
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The second term in the bracket is of the form 


_ _ — _ _ 

- cr 4 )r'^ 3 a;(r, - (74)0; —, ag - p 

V^o 

Using (3.20), this term can be distributed in /ii, h 2 , as terms of order O(r^). A similar 
argument holds for the third term. Indeed, we introduce a compensator 



■ a;(g,a)-a;(g,/3) 




, a ^ 13, 
a = I3, 


^i^,0!,(3) = ^a,/3(0 = 

allowing to rewrite this term as 

?]4(iy)r'^^('a3 -a4)Q I —,^3 - P,^4 - p] ■ 

Vo ) 

Again, using (3.20), this term can be distributed in /ii, /i2, as terms of order 0(r'^). 


(3.24) 


This allows writing the displacement map as a sum of four terms 

U(r,p) = -£o(M)( 1 + /io(r,p)) - C' 3 (M)£i(M)r'^^(l + /ii(r,p)) 

+ */l3r°'V(l + h 2 {r,p)) + K{My^ffu ( —,0^3 -p 

Vo 

with ho, hi of order OV). Moreover, K{M) = p 4 {y) — ri 3 {y){l — ei{M)) = Op{Mc). For 
p > 2, we conclude that the cyclicity is at most 3 by Theorem 5.12. 

For p = 1, we will prove in Theorem 5.13 that the cyclicity is at most 2. To this end, we will 
use that ^74(0) = —^3(0) = JI 3 and then that K{M) = + 0{iy)0p{Mc), in order to rewrite 

Uas: 

V{r,p) = -eo{M){l + ho{r, p)) - C 3 {M)ei{My^{l + hi{r, p)) 

+ *p3r'^^p{l + h2{r, p)) + *p3r'^^puj —, ag - p 

V^o 

with ho, hi and /ig of order OV). □ 


(1 + hoir, p)). 


(3.26) 



3.3 Finite cyclicity of the boundary limit periodic sets of 

Theorem 3.6. The boundary limit periodic set of {Hlfj (see Figures 1(c) and 5) has finite 
cyclicity. 


Proof. The proof is very similar to that of Theorem 3.3. The graphic occurs in the family (3.1) 
for B < I , which corresponds to I < a < 1, but we prefer to use the following equivalent un¬ 
folding inside quadratic systems (only parameters’ names are changed so that they play similar 
role as in Theorem 3.3) 

f X = -y + Bx"^ + p2y‘^ + hbX 27) 

\ y =x + xy + /i3p2 + p^y. 
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Figure 5: The boundary graphie through Pi and P 2 and the four seetions Sj and 11^,^ = 1, 2, in 
the normalizing eoordinates. 


The point P 4 (resp. P 3 ) is replaeed by Pi (resp. P 2 ). The quantity cxj is now given by cxj = 

The main differenee with Theorem 3.3 is that the transition from 112 to IIi is replaeed by the 
eomposition T~^ o D~^ oT o Di o Ti. The transitions Ti and are along the equator of the 
Poineare sphere and henee preserve the eonneetion (no translation terms). The saddle points Pi 
and Pr have inverse hyperbolieity ratios: = l/r,, = < 1. Henee, it is better to eonsider 

a displaeement map 

1/ : S 2 — )■ n^, V = T o Di o Ti o D 2 — Dr o Tr o Di o S. (3.28) 

The eomputation of S is the same as before. 

Computation of Ti and Tr. Tr and Ti are regular C'^-transitions with no translation terms. 
They can be computed in the coordinates (n, w) = (~yj y)- The transformed system in these 
coordinates is given in (6.1). The transitions take place along w = 0. Along this line, div = 
(3 — 2B)v — 2 /i 3 . Hence T)!(0) — T^'(O) = 0 (/r 3 ) = This property is preserved in the 

normalizing coordinates. 

Computation of T. The transition T in studied in (3.2). The line ?/ = —1 is invariant under 
= Hi. Hence, the constant term is of the form 

r(0) = eo(Tf) = *(/i4 - i^/ig). (3.29) 

Under the condition Eq = 0, we have div|y=_i = (2P + l)a: + /is — vJIs- Hence, 

r'(0) = ei{M) = */i 5 + 0(/i4) + 0(z/)0(/i3). (3.30) 

The equations (3.29) and (3.30) remain valid in the normalizing coordinates, and we call the 
corresponding coefficients Eq and £i. 
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The Dulac maps Dg, and Dr. We first localize the system (3.27) using coordinates (m, z) = 
(^, ^). The normalizing coordinates are of the form (m,, z), i G {f, r}. Then, 



(3.31) 


with C, a C^-function on monomials (see Appendix II). 

The Dulac maps Di and D 2 . They are given in Theorem 2.2. Since the connection along the 
equator is fixed, then the coefficient rji vanishes identically when cxo G N. 

Hence, the displacement map V (r, p) has the form 


V(r,p) = £ 0(1 + ^o(ap)) + + hi(r,p)) - + h 2 (r,p)). (3.32) 


This equation contains no resonant monomials since <72 + 7 ^ 1 as soon as B ^ 

We conclude that the cyclicity is at most two by Theorem 5.8. 


□ 


3.4 Finite cyclicity of ( 1 / 4 ) 

We now prove Theorem 1.2, i.e. that the graphic (/I 4 ) has finite cyclicity inside quadratic 
systems (see Figure 1(a)). 

Proof of Theorem 1.2. Such a graphic occurs for system (3.1) when Bq > 1, and its deformation 
in quadratic systems is given in (3.2). As usual, we should normally consider all limit periodic 
sets of Table 2. It was shown in [ 8 ] that a graphic through a nilpotent saddle point has finite 
cyclicity inside any -unfolding under the generic conditions that the return map P along the 
graphic has a derivative different from one and that the nilpotent saddle point has codimension 
3. But the only limit periodic sets of Table 2 for which we use the genericity hypotheses are the 
boundary limit periodic sets which have been treated in Theorem 3.3, and the intermediate and 
lower limit periodic sets of Sxhhl and Sxhh5. 

For these limit periodic sets, we only have Dulac maps of the first type as in Theorem 2.2. 
Hence, we can work with a 1-dimensional displacement map, which we take as 1/ : S 3 —)■ 114 , 
V = D 4 o S — T o D 3 (see figure 6 ). As before the sections Sj and Hj are parameterized by 
the normalizing coordinate Xi near Pj, which are chosen so that S and T are the identity in the 
center case. 

The technique is to write V in the form of a finite sum 



(3.33) 


for some cr > 0. The parameters are the same as in (3.9) and (3.10). We write little details since 
they are very similar to [7] . 

The intermediate graphics. For these graphics, the map V{x 3 ,p) is in x^. Under the 
condition /r 2 = /^s = 0 for a nilpotent saddle, (3.2) has an invariant parabola for = 0, which 
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Figure 6 : Intermediate and lower limit periodic sets of Sxhhl and Sxhh5: the four sections Sj 
and Ilj, i = 3,4, in the normalizing coordinates near P 3 and P 4 . 


is the only possible connection at a nilpotent saddle. Hence, T has a constant term of the form 
+ 0(/i3) + The constant term of the transition S has the form since JI 2 

respects the symmetry, and hence does not contribute to the breaking of the connection. 

When (Jo ^ N, this yields that the constant term £0 in the displacement map has the form 

Eq = *^4 + 0 {iy) 0 {^^) + 

When (To = P € there are additional terms 


z/ 


^0 


0'3 


Wo — W —,03 


Z/ 




z/ 




(74 


Wo ~ a; —,04 


t'o 


+ ^?4(a3 - a4)Po ( — ) u}{—,a2,-a4]u}{—,az 

\^qJ V^o / V^o 

+ ^4(03 - a 4 )Po ^ 


(3.34) 


—, 0^3, 0^4 

i^O 

In this expression rj^ — ri 4 = Op{Mc) and as — 04 = 0p(Mc)0(iy). Hence, in all cases we 
have 

£0 = *M4 + 0 (z/) 0 (/Z 3 ) + M 50 (M) + 0(iy)0p(Mc). (3.35) 


The linear term has the form iy°'^T'(0)— iy‘^'‘S'(0)). Moreover, S''(0)|p=o = 1 precisely when 
Jls = 0. Also, Lemma 6.2 shows that T'(0) — 1 = */i 5 + 0 (/i 4 ) + 0{^s)- Considering that 
0^3 — 0^4 = O^u), then 


z/'^4 ^ + (0-3 - (j4)a;(z/,(73 - (J4) = + 0 { v )). 
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This yields 


(*/^5 + 0(/i4) + 0(l^)0(/i3)) . 


(3.36) 


Now, because of the funneling effect, any nonlinearity on the side of T has a high coefficient 
in V which damps it. Hence, the only significant nonlinearities are on the side of S. We are sure 
that S is nonlinear when /I 3 7 ^ 0. This comes from the fact that the graphic belongs to a 
family of graphics. In the case of Sxhhl, this family ends in a lower graphic with a saddle 
point and its hyperbolicity ratio r is different from 1 precisely when /Ig 7 ^ 0 , yielding that 
S{x 3 ) = Co + C 1 T 3 + o{xl), with Cl 7 ^ 0, for graphics near the saddle point, and hence that S 
is nonlinear on the whole section S 3 . Then, for any graphic occuring for a value x^^, there exists 
n such that 7 ^ 0. Hence, C(")(T 3 ,o) = [cn, 37^3 + 0{v)Op{Mc)] = Sn- 

Moreover, for all graphics except a few isolated ones we have that n = 2. The same argument 
can be applied for Sxhh5 since the connection is fixed between the two saddles and the product 
of their hyperbolicity ratios is different from 1 precisely when /xg 7 ^ 0. Hence, we have written 
V under the form (3.33) with hi{x 3 ) = ^^(l + 0 (a; 3 )). We conclude to finite cyclicity by means 
of Theorem 5.8. 

The lower graphic of Sxhhl. The study is very similar and divided in two cases. When 
/Z 3 7 ^ 0, it was already shown in [ 8 ] that the lower graphic of Sxhhl has finite cyclicity. This 
comes from the fact that the hyperbolicity ratio r at the saddle point is non equal to 1 precisely 
when /I 3 7 ^ 0 , in which case we conclude to finite cyclicity because of the nonlinearity of S. 
Hence the difficult case is the neighborhood of/Z 3 = 0 since, for this value, tq = 1. In that case 
we reparameterize the section Eg by means of xg = Tg — co(M), so that xg = 0 corresponds to 
the unstable manifold of the saddle point on the blow-up sphere. Then, as before, we write V 
as a sum of terms: 


V (xg, M) = Eo^oixs, M) + fi^X 3 Uj{x 3 , T - l)/l3(x3, M) + £4X3/11 (xg, M), (3.37) 

with hi{ 0 , 0) 7 ^ 0. We conclude to finite cyclicity by means of Theorem 5.8. 

The lower graphic of SxhhS. Such a graphic occurs for /Z 2 > 0. Because the connection 
is fixed between the two saddles, the map S can easily be computed and has the form cq + 
CiXg + o(x 3 ), where r = 1 ^=3— is the product of the two hyperbolicity ratios. Again, we 

reparameterize the section Eg by means of xg = xg — co(M), so that xg = 0 corresponds to the 
unstable manifold of the right saddle point on the blow-up sphere. This allows writing the map 
V in the form 

fv"(x3) = eix\hi{x3, M) + JI^xlhr{x3, M), tq ^ N, 

[(/(xg) = Y 2 = 0 ^i^\hi{x 3 ,M) -h/Z3X3“w(x3,r-ro)/lr(53,Af), Tq G N, 

with /ij(0, 0) 7 ^ 0. We conclude to finite cyclicity by means of Theorem 5.8. □ 
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4 Appendix I — Hyperbolic fixed points 


We will consider germs of smooth family of 3-dimensional vector fields o- at (0) G with 
coordinates (m, n, y), which are quasi-linear of the form: 

{ ii = u, 

i) = -V, (4.1) 

y = -(^y + F^{u,v,y), 

where a is a parameter in a neighborhood of ao G M’*', and y a parameter in a neighborhood of 
yo in some Euclidean space. Moreover, = 0{\{u,v,y)\‘^) at the origin, for any value of the 
parameter {y, a). The system has the first integral: a = uv. 


4.1 Normal form 

It is possible to find local normal form coordinates for by a coordinate change preserving 
the coordinates u and v. More precisely, we have the following normal form result: 

Theorem 4.1. There exists a normalizing operator fif defined on each pair (X^ c,, k), where 
X^ o^ is a family as above and k G N*, such that, 

F) ifiki F{kfi G , 


where 

{u,v,y) {u,v,Y = 

is a parameter-depending change of coordinates of class defined defined for \a — ao\ < S, 
\y — yo\ < Ek, and \ {u,v,y)\ < pk, such that 0, 0) = Id, which brings to the 

following polynomial normal form of degree K{k): 

1. If CIO ^ Q : 

{ ii = u, 

V = -V, (4.2) 

y =-W + <Pf,A’'))r 

2. Ifao = I G Q, with (p, q) = I when q 1 : 

{ ii = u, 

V = -V, (4.3) 

with = 0 when ao ^ N ( q 1). 

The functions Tfii,u are polynomials of degree < K{k), with C°° coefficients in {y, a) 

and 0) = 0. 
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Proof. The proof is standard in the literature, and we only reeall the main steps. 

The degree K{k) can be determined algorithmically from the eigenvalues {1, —1, —cxo}. 
The number 6k is chosen sufficiently small so as not to introduce any new resonant terms of 
degree < K{k) for some a e [uo — 6k, <Jq + 6k]. 

The first step is to bring the system to normal form up to degree K(k) 

{ ii = u, 

V = -V, (4.4) 

z = P{a, /r, u, V, z) + R{a, /i, u, v, z). 

where P{a, /i, u, v, z) is a polynomial mu,v,zof degree K(k) containing only resonant terms, 
and R{a, /i, u, v, z) = o(|( m, v, This can be done by means of a polynomial change of 

coordinate 

K{k) 

y = z+ ^ aijiPffz^. 

i+j+£=2 
i-j+o-Q (£-1)^0 

Because this change of coordinate is tangent to the identity and contains no resonant monomial, 
then it is uniquely determined. 

The second step is to kill the remainder R in (4.4). For this purpose, we decompose R as 
R= Ri +with Ri = and i?2 = Each part is killed by the 

homotopy method. The details are exactly the same as in [5]. Again, this step is algorithmic. 

□ 


4.2 Properties of compensators 

This section is devoted to properties of different fonctions useful for the expression of the re¬ 
sults, and in particular the so-called compensators 0Ja{0 VLa,i 3 {i) defined in (2.9) and (3.24). 

First, we introduce the analytic function 


K{ri) 


1 , 77 = 0 . 


(4.5) 


The following Lemma gives some useful properties of k: 

Lemma 4.2. The function k is an entire analytic real function whose series is given by nff) = 
(thT)!- Rfollows that ^{rj) < K{rj) < for rj > 0. Moreover, K{ri) > 0, ^{rj) > 0, and 
> OJorallr] G M. 

Proof We have that n{r^) = 4(Er S - 1) = and then: ^(77) = ^y 

The inequalities ^(^7) < ^(77) < for 77 > 0, follow trivially. 

Clearly, ^(77) f 0 for all 77 G {0} and as k(0) = 1, it follows that ^(77) > 0 for all 77 G M. 
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Next, as ^{rj) = , any root rj ^ 0 of ^{rj) = 0 verifies that Comparing 

the series of these two funetions, we see that for rj g]0, 1[. The inequality ^{r]) > 0 

is trivially verified when r] > 1 . Finally, when 77 < 0, we put r] = —5, with 5 G The trivial 
inequality: > 1 + 5, for 6 G implies that for rj < 0. As ^(0) = we have that 

^{r]) >0 for all r] eR. 

To finish, since ^ij}) = any root 77 7 ^ 0 of ^(r;) = 0 verifies that = 

2 , or equivalently e~^ = 1 — 77 + \rf. Let g{rj) = e“^ — 1 + 77 — Let us show that 

F 2 V 

div) 7^ 0 foi" ^ 7^ 0 - Indeed, g'{r]) = —e~^ + 1 — 77 < 0 . The numerator of is —e^g'{g) 
and is positive for 77 7^ 0 . Hence, g'{g) < 0 for 77 7^ 0 , and since 5f(0) = 0 , then gg{g) < 0 for 
77 7^ 0 . As ^( 0 ) = we have that ^(77) >0 for all 77 G M. □ 

The following lemma gives the relation of ui defined in ( 2 . 9 ) with k, and interesting proper¬ 
ties which can be easily deduced using this relation: 

Lemma 4.3. We have that a) = —^{—a In,^) In,^. The compensator oj verifies the follow¬ 
ing estimates 

1. uj{^, o) < — In^ if a < 0 and a;(^, a) < —In^ if a > 0, and then 

a;(e,«) = 0(r'“'|lnel)- (4.6) 


2 . 

a) —)■ +cxd when {^,a) —>■ (0,0). (4.7) 

Proof Using properties of n given in Lemma 4.2, it follows that: 

1. If a > 0, i.e —aln^ > 0, then a;(^, a) = —/s:(—aln^) In^ is less than —In^ = 

-r“ine. 

2 . If a < 0, i.e —aln^ < 0, thena;(^,a) = —k(— aln^) In^ < —In^ (indeed, n is 
increasing, 7c(0) = 1 , yielding K{g) < 1 when 77 < 0). 

The estimate (4.6) follows from these two inequalities. In order to prove (4.7), we take any 
A > 0. 


1. If —aln,^ > —K, we have that /?(—aln,^) > k{—K), as n is increasing, and then 
> -k{-K) In^. 


2. If —a In,^ < —K (in particular a < 0), we have that 


a;(^,a) 


1 - ^ 1 - e"^ 

lal ~ lal 


from which (4.7) follows. 


□ 
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In parallel with the eompensator 11 introdueed in (3.24), we introduee the symmetrie fune- 
tion 




k{vi)-k( 5) r, =L S 

r^-S y T 

t^V)y V = S. 


(4.8) 


This yields 


^(C,a,/3) = /C(-aln^,-/31n0 In^^. 


The useful properties of f2(^, a, /3) are given by the following lemma: 


Lemma 4.4. lla,/3(0 = 0{^ '^In^.^), where 7 = max{|Q;|, |/)|}, and lla,/3(0 + 00 , when 

(C,a,/3)^(0,0’,0). 


Proof. To prove the two elaims, we just have to use the Mean Value Theorem for the function 
/C: there exists 6 e [ 7 ,5], such that /C( 7 ,5) = 

Let us begin by the first claim. Let us start with the case a > jd. Then/C(—aln^, —/)ln^) = 
^( 6 *), for some 6 G [—fd In^, —a In^]. As ^( 7 ) is an increasing function (see Lemma 4.2), we 

have that/C(-aln^,-/)In0 < If a < 0, we use that g(-alnO < f (0) = I 

to obtain that lla,/ 3(0 < ^ If a > 0, again using Lemma 4.2, we have that g(-alnO < 
g-aing _ ^(^) < ^ summarizc the two possibilities by 

writing that ^a,fsiO ^ soon as a > (d and ^ and |a| sufficiently small. Using 

the symmetry of ^a,i3iO can permute a and fd in the above argument to obtain finally that 
= ^(^-^In^O, where 7 = max{|a|, |/3|}. 

We now prove the second claim. By symmetry on a and fd it suffices to prove the claim for 
a > fd. As above, we can write that ^a,i3{0 = some 6 G [—/31n^, —aln^]. 

Now, we want to bound from below. Since ^ is increasing, 11 q,,/ 3 (.^) > ^{—fd In.^) In^ 

If /3 > 0, we just use that ^{—fd\n^) > ^(0) = to obtain that lla,/3(0 ^ ^ 

/) < 0 , we have to compute |(-;31na = liif) = g(lnfW). As 

we have that |(-/3hiO = and then: > (M'vS-y'-'vi , yieui„g ,hat 

^q:,/ 3(0 + 00 - This yields the conclusion. □ 


4.3 Transition along the trajectories 

Bfe want to study transition maps for (4.1), in the region Q = {u > 0, u > 0,} C near 
the origin. More precisely, let VL be a neighborhood of the origin in M^, and 11 C {m = mq}, 
for uq > 0, be a section. The neighborhood W can be chosen sufficiently small so that the 
trajectory starting at any point in VL fl {u > 0} reaches If for a finite positive time (in particular, 
lU n n = 0). We consider the transition o- from the points in lU n {m > 0} to the section 11. 

We will compute T^,cr, in the -coordinates given by Theorem 4.1. In this system of coor¬ 
dinates the family is the smooth family of polynomial vector fields (this means polynomial 
in (u, u, y) with smooth coefficients in (/r, a)). 
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We take 11 = [—^ 0 ,^ 0 ] x [ 0 ,^ 0 ] x {mq} for some Yq > 0,vo > 0. On 11, we replace 
the coordinate v by u = uqv, with u G [0, z/q = wot^o]- Then, we can write T^^„{u,v,Y) = 
{Y^^„{u,v,Y),v = uv). 

The expression of the Y -component Y^^„ is given by the following Theorem: 

Theorem 4.5. Let a = cf (a, v) = a + cind a = a(a, u) = d{a, u) — cxo, where is 

the polynomial family introduced in Theorem 4.1. The Y-component of the transition map 
has the following expression onW P\ {u > 0}.' 

1. Ifao ^ Q : 

Y^,,{u,v,Y)= {-) Y. (4.9) 

\UoJ 

2. Ifao = I e Q with (p, q) = 1, when (Tq ^ N : 


Yf,^^{u,v,Y) = 4-n)j, (4.10) 

where is the same as in (4.3) (in particular, = 0 when ao ^ N)- 

The function family in (4.10) is of order O and, for any 

integer I > 2, is of class C^~‘^ in (Y, , v, p, a). 

Proof The time to go from a point (m, n, F) G IF fl {m > 0} to the section If along the flow of 
is equal to — In Expression (4.9) follows trivially from the integration of the third line 
of the system (4.2). 

Then, from now on, we will assume that ao G Q and we will study the integration of the 
system (4.3). The trajectory through the point (m, v, Y) is equal to (we*, ve~^, Y(t)) where F (t) 
is solution of the 1-dimensional non-autonomous differential equation: 

Y(t) = -dY(t) + <h^,.(z/,M"eP*F(f)'?)F(f) + e-^Vp^,.(z/), (4.11) 

with initial condition F(0) = F. 

In order to eliminate the linear term in (4.11) we look for F (f) in the form F (f) = e~'^*'Z(t). 
As Y(t) = e“'^*F(t) — dY(t), and letting d = ^ 4- a, 'we obtain the following differential 
equation for Z(t) : 

Z = e-''“VZ‘')Z + e“Vp^,.(z/), (4.12) 

with initial condition Z(0) = F. Note that the term in o- is only present when q = 1. 

The 1-dimensional non-autonomous differential equation (4.12) is smooth in (t, Z, a, u, u, v, p) 
and can be integrated for any time f G [0, — In ^]. If Z(f) is the solution of (4.12) with initial 
condition Z(0) = F, we will have that 

Y^^,(u,v,Y)= (-)"z(-ln-). (4.13) 

\Uq/ V Uq/ 
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The above expression has to be studied for m > 0 (we extend V along [u = 0} by 
f, Y) =0). We first study the integration of (4.12). 

To begin, it is easy to get rid of the term in (4.12). Let us consider the analytic 

function 


0(f, a) 


—1 / r\ 

-, a ^ U, 

t, a = 0. 


which verifies 0 = e“*. We have that 0(f, a) = tK,{at) and then uj{^, a) = 0(— In^, a). 

Putting Z{t) = a) + Z(f), we see that Z{t) is the solution of the differential 

equation 


Z = + Z)-?) (r;%,.(i/)0(f, «) + Z), (4.14) 

with initial condition Z(0) = Y. 

As 0) = 0, we can write where is a smooth function. 

Now, let us notice that e"* = 0 = 1 + aQ. Moreover the map t —>■ 0(f, a) is invertible (for any 
a). Then, we can change the time t by the time 0 in the differential equation (4.14). We obtain 
the new equation 

d/Z 

= u^H^Q, Z,u,v,h',a, fi,a) (4.15) 

with 

H={1 + aQ)-^^+'^\vPr]Q + Zy+^H^^„(u,uP{l + aQ)-^{vP7]Q + Z)'?), (4.16) 

where r] = Let fh (0, Y, u, v, z/, a, /r, a) be the solution of (4.15), with the “time” 0. 

Up to now, 0 is seen as an independent variable; in particular it is independent from a. For 
^ ® yielding 


Z 



= ^ 



, y, M, V, o, a, p, a 



(4.17) 


and then, the computation of v, Y) reduces to the computation of 4/ aj, Y, u, v, p, a 

One difficulty in the study of 4/ j, Y, u, v, u, a, p, cxj is that u aj — )■ +cxd if 

u —0. To overcome this difficulty we will exploit the fact that the right hand side of (4.15) is 
divisible by vP. 

We first study the differential equation (4.15). We put u = and change the time 0 by 
the time t = UQ (and not just by uQ, as it could seem more natural). The equation (4.15) is 
replaced by the following equation 


d^ 

dr 




(4.18) 
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where H is given by (4.16). Let G be the right hand side of (4.18). It is smooth for U > 0, 
but since it is function of ttu, it is not well-defined in a whole neighborhood of the point 
{(r, Z, f/, V, V, a, /r, a) = (0, 0, 0, 0, 0, 0, /tq, o'o)}- Fortunately, we only need to integrate (4.18) 
in a closed domain V\ 

Definition of V. The domain V is defined in the space (r, U, Z, v, u, a, /r, a) defined by 

1. f7 e [0, t/i], |a| < tto and r G [0,Uuj{^,a)], where Ui,ao > 0 are chosen arbitrarily 
small (the time r = Uoj( — , a) corresponds to the time f = — In — = — / In —), 

2. (Z, V, u, a, yU, a) G A, an arbitrarily small closed neighborhood of the value (0,0, 0, 0, yUo, cto)- 

We want to prove that G is of class on V. We will first prove a technical lemma about the 
partial derivatives of the function G. Let us denote by dmG any partial derivative of G corre¬ 
sponding to a multi-index m = (mi,..., m^) associated to the variables r, U, Z, v, z/, a, yU, a 
and the coordinates of yU. Let \m\ = mi + ■■■+ iris he the degree of m. We will note by 5, a 
strictly positive number, which can be made arbitrarily small by appropriately choosing Ui and 
A. We have the following: 

Lemma 4.6. Let be ^ as above. Let m be any multi-index such that \m\ <1 — 2. Then, 
for any 5 > 0, there exists a domain V as above, such that on the restriction to the domain V 
we have that 

= 0(f/P'-H-i-^). (4.19) 

Proof. Recall that G = U^^H, where H is given by (4.16) and 0 is replaced by The proof is 
straightforward, but rather tedious, and we just give the main steps. First, let us notice that on 
V we have that, for any s G Z: 

(l + a^y = (1 + aoy = = 0([/-l"^"l). (4.20) 

Also, using Lemma 4.3, we have that: 

1 = 0 = K(at)t < lnU\. 

/ \ -(9-1-1) _, 

These estimations imply that + and ^ have an order 0(U ). As if is bounded 

on V, we have that G = This is the expected result for m = 0. 

Next, we use the expression of the partial derivatives of G, in terms of the functions 0, 

(1 -f 0:0)“'^ or (1 -f q; 0)“('?+^) and the partial derivatives of if^^cr, evaluated on V (these partial 
derivatives are bounded on V). We have for instance that: 

^(1 + a0)-'? = -qlail + a0)-('?)^ = OifJ-^-^). 

As (1-|-Q;0)“'^ = 0{U~^), we remark that the order in U has discreased by one unit (modulo 
an order in 5). 
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It is easy to see that this observation can be generalized for any partial derivative: the 
previous order in U decreases by one unity for each first order partial derivation (modulo an 
order in 5). 

Then, starting with G = for m = 0, the estimation (4.18) for any multi-index 

m follows direetly by reeurenee from this fall of order (let us notiee that, in a symbolie way, we 
have: “5 + 5 = 5”). □ 

End of the proof of Theorem 4.5 

Lemma 4.6 says that eaeh partial derivative dmG ean be extended eontinuously on r = 
(/ = 0 by giving it the value zero at these points. Then, as the funetion G is smooth on 
'D\ {t = U = 0}, the restrietion of (5 to P is a funetion of differentiability elass on the 
whole domain "D, including the points on {t = U = 0}, when we give to eaeh partial derivative 
of G or order less than I — 2 the value 0 at these points. Let B he a elosed neighborhood of 
(0, 0, 0) in the (r, a, (7)-plane, eontaining the elosed set 

{(r, a,U) \ t e [0, —/[/In ^], |a| < a^, U E [0, Ui]} 

Uo 

that we have introdueed above in the definition of V. The elosed domain V is eontained in the 
neighborhood Ax B. Using the Whitney Theorem for the extention of differentiable funetions 
(see [M] for instanee), we ean find a C^“^-funetion G on a +l x sueh that G\^ = G (here, this 
extention ean also be easily eonstrueted by hand, in an elementary way). 

For times r G [0, —lU In the flow \I/(r, Z, U, v, z/, a, /i, a) of the differential equation 
(4.15): = G eoineides with the flow Zp,a) of the differential equation 

^ = G. This equation is of differentiability elass on Ax B^ as well as its flow 'F. 

In partieular, we have that 

z(^ — \n —j —’ '^) ’ d) ) (4.21) 

is a C^“^-funetion of (F, U, , v, u, a, /r, a), i.e. is a C^“^-funotion in the variables 

(y, ut ,uTu}{f^,a),v,u,a, ^, 0 ^, a function which is defined on a neighborhood of the point 
(0, 0, 0, 0,0, /io, (To). We ean replaee a (outside uj) by its expression in (cr, u) and u by uv to 
obtain finally that z(^ — \n ^ j is a -funetion of ufiuTp, As Z{0) = Y, 
we ean write 

z(-ln-) = Y + (j)^AY,u,v), (4.22) 

V Uq/ 

where ^ 

^[/a;(^^,aj,F, f/,n,z/, a,/r, (T^ - F (4.23) 

is a C^“^-funetion of {Y,ufiuTuj{^,Q.),v, . Finally, eolleeting the different terms in 
(4.13), (4.17),J4.22) and (4.23), we obtain the expression (4.10) in Theorem 4.5, for the transi¬ 
tion funetion v,Y). 
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We can estimate from the differential equation (4.14) for Z(t). If G{t, Z, u, v, u, a, a, jj) 
is the right hand side of (4.14), we have that G = 0{u ^on the domain V defined 

above. As t < — In;^ on P, then 0(^,0;) < yielding G = 

From this estimate of the order of G, it follows that 

0^,^ = Z^-ln—) - Y = ,a)| Inwl), 

V Uq/ \Uq / 

which is the estimation in the statement of Theorem 4.5. □ 

4.4 Transitions between sections 

Theorem 4.5 gives the expression of the transition = (z/, F^ o-), starting from any point 
(m, n, Y) in the domain hF fl {m > 0} and landing on a section 11 C {m = mq}, for some mq > 0 
(we can extend trivially to the whole neighborhood W by taking Y^^o-iu, n, 0) = 0). We 
apply this to get Theorems 2.2 and 2.3 after changing (m, v) i— )■ (r, p). 

Discussion of Theorems 2.2 and 2.3. A previous version of Theorems 2.2 and 2.3 was given in 
Theorems 4.10 and 4.14 of [8]. It is interesting to compare their proofs and formulations with 
the proofs and formulations in the present paper. 

1. The proof in the present version is unified: Theorem 4.5 gives a formula for a global 
transition from any point in a 3-dimensional neighborhood W, formula which is easy to 
restrict on the two different types of section E. Next, the proof of Theorem 4.5, even if it 
is based on the same normal form, is much shorter than the proofs of Theorems 4.10 and 
4.14 given in [8] . The reason seems to be that in [8] the transition function Y and its partial 
derivatives are directly estimated by a variational method. In the present paper, we have 
replaced the 1-dimensional non-autonomous differential system: Z = G, which is not 
defined in a neighborhood of the point {(r, Z, U, v, z/, a, p, a) = (0, 0, 0, 0, 0, 0, /xq, 

by a differential equation: Z = G, differentiable on a neighborhood of this point. As 
a consequence, we obtain almost without computation that the function 0^ g. is differen¬ 
tiable (in terms of fractional power and a compensator of some variable). In fact, the 
heavy computations made in [8] are replaced by an implicit use of the Cauchy Theorem 
for differential equations. 

2. We can compare the statements in [8] and in the present paper. We restrict the comparison 
to the only non-trivial case: do G Q. The transition function called here is given by 
the formula (4.11) of Theorem 4.10 of [8]. We can observe that it is quite similar to 
the above formula (4.20), up to the changes of notations. The same remarks are valid 
for the transition of type II which is treated in Theorem 4.14 in [8]. The only important 
difference is in the form and properties of the function which is called (j) or 9 in [8]. 
We will comment on this in the next items. 

3. The function 0^ in Theorem 2.2 is of order 0(z/P+'?"a;'^+^ aj | In z/|). This order has 
to be compared with the order given for the function 0 in Theorem 4.10 of [8] which is 
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exactly the same order for a < 0, but equal to aj | \nv\) for a > 0. This 

minor difference is probably due to the difference in the method of proof. It is less easy 
to compare the order of in Theorem 2.3 with the order of 9 in Theorem 4.14 of [ 8 ]. 

4. In Theorem 4.10 of [ 8 ], 0 is a -function of and other variables. Since u —>■ 

+00 for z/ —)■ 0, this means that the domain of 0 has to be unbounded. This implies that 
it is not possible to deduce directly the order of the partial derivatives of 0. This order 
is obtained by using variational methods and heavy computations. On the contrary, the 
formulation given in Theorems 2.2 and 2.3, permits a direct deduction of the order of any 
partial derivative of 0^ ct- Let us show this on an example for a transition map of type I. 
Considering any I G N and observing that 0^,^ is of order we can write 

0 /i,O' *0/^,t7) 

where is a C*“P“^-function in (Y, z/^/0 aj, /r, a). 

As a consequence any partial derivative of in terms of Y, /i, a, of degree less than 
/ — p — 3, is of order ^). Taking into account that we can take I arbirarily large, this 

order in very similar to the order obtained in Theorem 4.10 of [ 8 ]. 

5 Appendix II—Counting the number of roots 

5.1 Differentiable functions on monomials 

We come back to the notations of Section 3: r, p are variables defined in a compact neighbor¬ 
hood A of (0, 0) in the first quadrant Q = {r > 0, p > 0}. We will always choose A to be a 
rectangle [0, ri] x [0, pi], in order to have connected curves /y = {(r, p) & A\rp = 0}. In the 
following definitions we will use also compensators oj^ and depending on other param¬ 
eters 7 , 8. We will often use the shortened notation ce^y, > ^ 7,5 • Moreover, 

changing r to we can of course suppose that tq = 1. 

We consider a multi-parameter A in a compact neighborhood of a value Aq in some eu¬ 
clidean space 8. The neighborhood B will be chosen sufficiently small to have the desired 
properties. 

We also consider functions which are differentiable on real powers of r, p and compensators 
in r. We give a precise definition of this notion. 

Definition 5.1. 1. A primary monomial (monomial in short), is an expression M = r“, p^, 

or u^{r)~^ where a,b,c,d,G cind 7 , 71,72 are smooth functions 
of X. Moreover a, 6, e are strictly positive and 7(Ao) = 7i(Ao) = 72 (^ 0 ) = 0 (we can 
have y = a or (3 and ( 71 , 72 ) = (a,0)). For instance, rs, are primary 

monomials but not r°‘ or uj^- 

A monomial M defines a X-family of functions M{r,p, X) on Q = {r > 0, p > 0}, M 
is smooth for r > 0 and, by Lemmas 4.3 and 4.4, it can be extended continuously along 
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{r = 0}); we have that M(0, 0, Aq) = 0 (i.e. M = o(l), in terms of some distance of 
(r,p, A) to (0,0, Ao)). 

2. We say that a function /(r, p, \) on Ax B is a C^-function on the monomials Mi, ..., Mi 
if there exists a C^-function /(^i,..., A) defined on Ax B, where A is a neighborhood 
o/O G such that f{r, p, A) = /(Mi,..., M;, A). If the number of monomials and their 
type is not specified, we just say that f is a -function on monomials. 

Clearly, the space of C^-functions on monomials, defined on ^ x is a ring. The classical 
theorems of differential calculus (Taylor formula, division theorem and so on) can be extended 
to these functions by applying them to the function /. Since the differentiability class k is finite, 
there will be falls of differentiability class in these operations: Lemma 5.3 is one example. For 
this reason, we will consider functions / with the property to be -functions on monomials, 
for any k G N (but with a choice of monomials and a size of the neighborhood A x B that 
may depend on k). The functions and introduced in the 

statements of Theorems 4.5 are, 2.2, and 2.3 are examples of C^-functions on monomials for 
any k, which use only the single compensator cJq. The functions hi entering in the expression 
of the displacement map V in Section 3 are using other compensators ui^, and also kla,i 3 - 


5.2 Procedure of division-derivation for functions with 2 variables 


Notation 5.2. In this section, h{r, p. A) = o(l) will mean that h(0, 0, Aq) = 0. 


We want to bound the number of roots of an equation {l/(r, p. A) = 0} along the curves 
lu = {rp = z/ I (r, p) G .4,}, for z/ > 0 and a neighborhood Ax B sufficiently small. The 
function V is expressed using C*^-functions on monomials. To obtain this bound, we will apply 
Rohe’s Theorem, and to this end we will use recurrently the Lie-derivative Lx of V hy the 
vector field 


y 9 d 

- %■ 


(5.1) 


Hence, we need some properties of acting on C^-functions on monomials. It is easy to see 
that: 


' Lxr^ = ar^, 

Lxp'" = -bp\ 

Lxio-i = — (1 + 

~ (^71 3 “ 72^71,72)- 


(5.2) 


From this, it follows that 


Lemma 5.3. If f is a C^-function on monomials, then Lxf is a ^-function on monomials 
and Lxf = o(l). 


Proof If M is any monomial, LxM is a linear combinaison of monomials. Then, Lxf — 
^LxMi, is a C^“^-function on monomials and, since each monomial is o(l), this function 
Lxfisa\soo{l). □ 
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For the procedure of division-derivation we will need more general monomials than the 
admissible ones: 

Definition 5.4. 1. A general monomial is an expression M = rij where i 

and j belong to finite sets of indices. The coefficients a, b, Ci, dj, as well as the internal 
parameters of the compensators Ui, Qj, are smooth functions of X (without any restriction 
on sign). Let a(Ao) = a°, b(Xo) = 6 °. 

2. A general monomial is resonant if (in this case the “polynomial” part r°‘°of 

M reduces to the first integral z/“ ). Seen as a function of (r, p, X), such a monomial is in 
general not defined for r = 0 and p = 0 . 

Remark 5.5. An interesting property is that if M is a general monomial, then M~^ is also a 
general monomial. 

Notation 5.6. For convenience, ifui = a;(r, 7 *) we will use the contracting expressions: u = 
7 = hi)i, c = Yl-ujf = uf Y.i'liCi = 1 C. 

A first easy result, which will be the principal tool in the proof of Theorem 5.8 below, is the 
following: 

Lemma 5.7. We consider an expression f = M(1 + h) where M = is a general non¬ 

resonant monomial without Vt-factor and h is a -function on monomials, of order o(l). Then, 
on a sufficiently small neighborhood B, we can write: 

Lxf = {a-b + yc)M(l + g), (5.3) 

with g, a -function on monomials, of order o(l). 

Proof. We have that Lxf — LxM(l -y h) + MLxh. Using the formula of derivation for oj, we 
obtain that LxM = (a — b -\- yc + cuj~^)M. As M is non-resonant, we have that a° — 6 ° 7 ^ 0 
and, if is a sufficiently small neighborhood of Aq, we will also have that a — 6 -f 7 c 7 ^ 0 on 
B. Then, we obtain that: 

Lxf = (a — b-\-'yc)(l-\ - - -'j M(\ -\- h) -\- MLxh. 

V a — b + yc) 

We can write this expression as Lxf = (a — b + yc)M(l -f g), with 

cuj ^(1 h) Lxh 

g = n-\ -;-. 

a — b -y yc 

It follows from Lemmas 5.3 and 5.7 that p is a -function on monomials, of order o(l). □ 

We want to use the algorithm of division-derivation in order to prove the following result: 
Theorem 5.8. Let V (r, p, A) be a function on A x B (1 {r > 0, p > 0}, of the form 

i 

V (r, p, A) = ^ Ai{X)Mi (^1 + gfr, p, A)), (5.4) 

i=l 

where: 
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1. the leading monomials Mi = are general monomials, without VL-factor (u = 

Ci = {cl)j with i G J, a finite set), 

2. the functions gi are C^-functions on monomials, with k > I, and of order o(l), 

3. the functions A* (A) are continuous, 

4. the monomials MjM~^ for i j are non-resonant, i.e. 

a° - a° - 6° + 6° 7 ^ 0 for i j. (5.5) 

Then, if Ax B is chosen sufficiently small, 

i) either the function V has at most I — 1 isolated roots counted with their multiplicity, on each 

curve ly = {rp = a} <Z A, 

ii) or V is identically zero. 

Proof We suppose that V is defined for A G (some neighborhood of Aq) and we define the 
following elosed subsets: 

S, = {A G I A(A) > AfiA),Vj = l,...,l}. 


Of eourse we have B = UiBi, and it is suffieient to prove the result for any Bi (and B 
suffieiently small). Then let us piek any i = 1,..., Z. By reordering the indiees, we ean suppose 
that we have pieked i = 1. 

The algorithm of division-derivation eonsists in the produetion of a sequenee of funetions: 
Vo = V,Vi,..., V^_i, sueh that eaeh Vj is a summation similar to V but only onl—j terms, and 
is defined on a smaller neighborhood A^ x B^ of (0,0, Aq). 

To define Vi, we first divide V hy Mi{l + gi) (a division step). This is made on a neighbor¬ 
hood A^ X B^ <Z Ax B chosen such that 1 -|- gfir, p, A) 7 ^ 0 for all (r, p, \) & A^ x B^. On this 
neigborhood we consider the function: 


V 

Miil + gf) 


k 

Ai + AiMiM.^ ^ -|- g^ , 

i=2 


where the function gi, defined by 1 -|- p* = is on monomials and of order o(l). 

Next we apply the operator Lx (a derivation step). Since the monomials MiMf^ are non 
resonant for i 7 ^ 1, we can apply Lemma 5.7 to obtain the following function Vi on A^ x B^: 


Vi = Lx 


V 

-Mfil + gi). 


i 

^(oi -ai-bi + bi)AiMiMf^(l + g]{y,z)'^, 

i=2 


with the function g}, on monomials and of order o(l). The effect of the derivation is to kill 
the first term Ai, thus reducing by one the number of terms in the summation. Except from this 
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fact, the terms of the summation are eompletely similar to the ones in V, but with the funetions 
Ai replaeedby (ai — ai — bi + bi)Ai, and the monomials Mj replaeed by the monomials 

For the recurrenee step of order j + l = — l,we assume that we have a function: 

I j 

m=l 

defined on some neighborhood x with functions gl, on monomials and of order 
o(l). As in the first step from V to Vi, we divide Vj by + (y'j+i j, whieh is possible 

on some neighborhood x B^'^^ C A^ x B^, and next apply the differential operator to 
produee a funetion 

I i+i 

2 = j +2 m=l 

where the gj^^ are on monomials and of order o(l). 

Performing the I — 1 steps of the reeurrence, we end up with a function 

Vi-i = {ai — bi — ai + bi) ■ ■ ■ (ai — bi — a/_i + + gj^, 

where gl is on monomials and of order o(l). 

As gl = o(l), and at least on monomials, we can choose a last neighborhood A x B^ C. 
X sueh that the function 1 + gl is nowhere zero on it. We restrict now \ E Wi = 
B^ n Bi- On this set we have the following alternative: A;(A) 7 ^ 0 or Ai(A) = ■ ■ ■ = Ai{\) = 0. 
In the last case, the funetion V is identical to 0 and has no isolated roots. 

Then we just have to look at values A where Ai(A) 7 ^ 0. For sueh a value of A, the funetion 
V;_i itself is nowhere zero on A x W. Consider now any eurve 1,^ in A. Reeall that the derivation 
L;!:' of a funetion G eorresponds to the derivation of G along the flow of X and that G is an orbit 
of this veetor field. Then, as Vi-i is equal to the derivation of V 1 - 2 , up to a non-zero funetion, 
Rolle’s Theorem applied to Vi_ 2 , implies that the restrietion of this function to G, has at most 
one root (let us notiee that G is eonnected!). The same argument based on Rolle’s Theorem can 
be applied by reeurrenee to obtain for each J<1, that the funetion has at most j — 1 roots, 
eounted with their multiplieity. Finally, the funetion V has at most I — 1 roots eounted with their 
multiplieity on c n A, for A e Wi. 

We obtain the result by considering in the same way the different subsets Bi. □ 

Remark 5.9. 1. Even if V is a summation on admissible monomials, it is clear that, in 

general, the division step may produce general monomials. This is the reason why we 
begin with general monomials in (5.4j. 

2. Using the first integral rp = u, we can rewrite the leading monomial Mj in the form 
Mi = We call Mi = a reduced monomial. The sum (5.4) may be 
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written in reduced form, with Pi = Oi — bp 

i 

V(r, p, A) = ^ iy^^Ai(A)rP^uj^^ (^1 + gi(r, p, A)), (5.6) 

i=l 

3. The non-resonance condition (5.5) in Theorem 5.8 is equivalent to the condition that the 
Pi(Ao) = Pi in (5.6) are two by two distinct. Up to a change of indices and a reordering, 
we can suppose in this case that p\ < p® ''' < P?- some of p^ may be 

negative, and also that one of them may be equal to zero. 


5.3 The results of finite cyclicity for the boundary limit periodic set 

We now want to apply Theorem 5.8 to the displaeement funetion V in the text. We write 
ds = (To + a. After putting this funetion in the redueed form (5.5), we have the following. 

1. In the case ao ^ N, the function V is given in (3.19) and we have the sequence of 

monomials: j.cro-i+aj allows applying Theorem 5.8, yielding that the 

boundary limit periodic set is at most 2. 

2. In the case ao = p E N, the function V is given in (3.25) or (3.26), and the sequence of 

monomials is: {1, r^cco,}. We have two resonant leading monomials when 

p f 1, and even 3 when p = 1. Theorem 5.8 does not apply in none of these cases. 


Hence, we give a direct proof for ao E N, using exactly the same procedure of derivation- 
division as in Theorem 5.8, but based on a more refined estimation than the formula (5.3) used 
to prove Theorem 5.8. Recall that the parameter was called M in this context. It will not be 
sufficient to consider the leading reduced monomials for M = Mq, and we will have to look 
more precisely at the form of certain remainders. 

We need the following result: 


Lemma 5.10. 




^a;„(^l + 0(/)) =+ 0(r^)), 


(5.7) 


Proof. We have that L 


A’ 


{1 + O(r^) 


= L 


A 


ty 

r OJa 




r°‘U}a 


r"r As r " = 1 -f acUa, we have that 


(1 + 0{r^)) + r°'UaO{r^). Now, 
= —r". Since 


r°‘UJa 


r°‘ijjaO{r^) is of order 0{r^) (for a smaller S), we obtain (5.7) by grouping the terms. 


□ 


Remark 5.11. The formula (5.7) is wrong in general if we replace the remainder by the more 
general remainder o(l). Let us consider for instance the expression f = r"a;a(l + p). We have 
that Lxf = —r"(l + p) — r'^UaP = —r“(l + p + ce^p). The term ojaP is not of order (^(1). 


Let A, B be neighborhoods defined as above. First we have the following result when ao f 1'. 


36 












Theorem 5.12. Consider the case ao = p E N, with p ^ 1. Then the cyclicity of the boundary 
limit periodic set is at most 3, namely for sufficiently small neighborhoods A and B, the equation 
V(r, p, M) = 0 has at most 3 roots, counted with their multiplicities, on each curve ly C A. 

Proof. Recall that the displacement map V is given by 

V(r, p) = *eo{l + ho) + + hf + *p 3 Z/rP-i+“(l + /is) + *K{M)Ar‘^Uo^. (5.8) 

The sequence of leading monomials in (5.8) does not verify the condition of non-resonance. To 
overcome this difficulty, we will use that there is no remainder in the last term, and that ho is of 
order O(r^). For hi and /is, it will be sufficient to know that they are o(l). 

As in the proof of Theorem 5.8, we define the partition B = BiU B 2 U in terms of 

the coefficients in (5.8). At each step we will have to restrict the size of B. We will not recall it. 

As the three last leading monomials in (5.8) are o(l), the cyclicity is trivially 0 when M G 
Bi- We suppose now that M E B2C B^U B4. Using (5.7), we obtain: 

= *eirP+^{l + gi) + + ps) + *K{M)Ar^. 

l + ho 

Now, the sequence of leading monomials r"} verifies the condition of non¬ 
resonance and we can apply Theorem 5.8 to Then, this function has at most 2 roots, 

and the function V itself has at most 3 roots, when M E B 2 U B 3 U B 4 . □ 

Finally, we have 

Theorem 5.13. Consider the case (Tq = 1. Then the cyclicity of the boundary limit periodic set 
is at most 2. 


Proof We can start with the formula (5.8) which is valid for any p eN. Moreover, for p = 1 we 
have that K{M) = p^fu) — p 3 (z/)(l -|- £ 1 ) = + 0{iy)0p{Mc). This is a direct consequence 

of the fact that the linear part of the system at the points P 3 and P 4 is given, up to a constant, 
by f = r, p = —p, y = —a{y -f psp). Then, we can split the last term in (5.8) as the sum 
*P 3 ur°‘Ua + vr°‘ojaO{v)Op{Mc). The second term gives contributions of order 0{r^) in ho, hi 
and / 12 , and produces a remainder /13 of order 0 {r^) for the last leading monomial r^Ua- 

Then, for p = 1, the displacement map V takes the form: 

V(r, p) = *£ 0(1 + ^ 0 ) + -f hf + *p 3 Z/r"(l -f / 12 ) -f *p 3 iyr'^Ua(l + /is) (5-9) 


The sequence of leading monomials in (5.9) does not verify the condition of non-resonance. To 
overcome this difficulty, we will use that ho and /13 are of order 0{r^). It will be sufficient to 
know that hi and /12 are o(l). 

As in the proof of Theorem 5.12, the cyclicity is 0 if |eo| > max{|ei|, IpsI}. 

Otherwise, let us consider Using (5.7), we have that 


U 




= *£ir^+"(l -f pi) -f *p3i/ ar"(l -f /12) -f r°Lxh2 + *p3Z^r-"(l -f go), 


l + ho 
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with g 3 of order 0{r^). Grouping the different terms, we obtain 


i/ 




_ ^ 1 +a 


* £i(l + gi) + */i3p(l + =t:a + g^) 


1 + ho 

where g 4 = *ah 2 + Lxh 2 + go is of order o(l). Now, the sequence of leading monomials {1, p} 
verifies the condition of non-resonance and we can apply Theorem 5.8 to This 

function has at most 1 root, yielding that V itself has at most 2 roots, if l^ol < max{|ei|, l/isi}. 

□ 


6 Appendix III 


Lemma 6.1. The parameter function Sq in the expression of the displacement map V has the 
form {3.16) for system (3.3). 

Proof Since the system has an invariant parabola for po = P 2 = Ps = P 4 = 0, it suffices to 
make the calculation for po = P 2 = Pa = Ps = 0. The system is integrable when p 4 = 0, with 
integrating factor (1 + y)^. Hence, it suffices to show that the following Melnikov integral is a 
nonzero multiple of p^. Indeed, 

□ 

Lemma 6.2. The parameter function ei in the expression of the displacement map V has the 
form {3.11) for both systems (3.2) and (3.3). 



Proof It has been proved in [3] (see for instance Theorem 3.5) that it suffices to show that 
J divdf = *p 5 along the invariant parabola when all parameters but ps vanish. Two different 
calculations are needed for the cases (3.2) and (3.3). In the first case, the invariant parabola is 
given by (3.4). Then, 


div dt = lim 


{2B + l)a: + (1 - B)po 


dx 


Xo^oo -y Bx'^ + BpoX 

1 + B{Xo + {B — l)p5)^ + ©(ps) 


= lim 

JVq —^00 


(2S + l)ln 


1 + B{Xo — {B — l)p5)^ + o(p5) 

_l_22^3/2(i^ — l)/i 5 ^arctan (^y/B{Xo + 0(p5)) j — arctan (^y/B{—Xo + 0{p^)) 
= 2B^h^{B - l)7rp5 + o(p5). 

The second case of (3.3) is easier since the invariant parabola y = ^ is independent of 

P 5 . Then 

dx 


div dt = 




2p5 


-f 1 


= 27rp5. 


□ 


38 











Lemma 6.3. The second derivative of the map S = p-F(0, p), where F is defined in (3.18) is a 
nonzero multiple ofjl^. 

Proof. We first localize the system (3.2) at the nilpotent point at infinity using the coordinates 
{v,w) = after mutiplication by w, this yields 

n = te + (1 — B)v^ — P 2 — dsv + vw{{3B — l)p^ + pf) + 

22 ( 6 - 1 ) 
w = vw — psw — (1 — 2B)p^w + vw . 

A similar localization can be done for (3.3). We now let the blow-up (u,w) = {rx,r‘^) for 
w > 0 , and we consider the restriction of the blow-up system to the (p, a;)-plane for r = 0 , 
(after multiplication by 2) 

P = -P{X-P3P) = P{P,X), 

X = 2 + {1 — 2B)x^ — 2jl2p^ ~P 3 xp = Q{p-,x). 


Note that this system is the same for (3.2) and (3.3). The singular points occur at T = ±/3 with 
fi = ^ 2^11 • We localize at P 3 using x^ = fi — x and at P 4 using X 4 ^ = fi + x. Hence, the 

system at P4 is obtained from that at P3 through (xs, fi) i-)- (—X 4 , —fi). The map is between 
two sections {xi = Xq} in the normal form coordinates x, near Pi and we take Xq small. 
The section {X 4 = Xq} (resp {X 3 = Xq}) has equation x = ffip) = —Xq + 0{p) (resp. 
X = / 3 (p) = xo -f- 0(p)). A formula for the second derivative was given in [ 8 ] (Proposition 
5.2), namely 


5"(0) = 5'(0) 


2 1 /I(0)S'(0) ( (0, /.(o)) - I'M (F I (0, / 3 ( 0 )) 


rhW /pii 

/ 

'/ 3 ( 0 ) V V 


P'O' 

!-^( 0 ,x) ) exp 


) ( 0 , x)dx ) dx 

'fm \ Q 


( 6 . 2 ) 


Here, 5"(0) = 1. We call the three terms in the bracket 2Ji, 2/2 and I 3 . Let us first consider J 3 . 


r, 1 r 8S-5 

l 3 = 4jlfi2 + {l-2B)xl)^(^-^B^ / (1 - Px^)(2 + (1 - 2P)x^)2(i-2-B)dx. (6.3) 

JXQ 

There are two different cases for /j(0) depending whether Pq = f or not. 

The case Bf) = |. In this case, the singular point has equal eigenvalues and a Jordan nor¬ 
mal form for nonzero fi^. Hence, the change of coordinate to normal form is tangent to the 
identity and /3(0),/4(0) = 0 {jl 3 ) 0 {Xo). Also the integral part of I 3 in (6.3) is equal to 
—2 ^|xo — In j 0. The result follows in that case. 

The case Po 7 ^ I* In this case, the change of coordinates to normal form is given by x = 

~ (^3 - 3 ^p) + C'(|(p,X 3 )n for P 3 (resp. x = -fi + (^x^P + 0 (|(p,X 4 )n for 

P 4 ), yielding /'(O) = ^(1 + 0(Xo)); 
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2/i + 2/2 — [[+]]/X 2 

As for the integral part in J3, it is given by 


Xq 


3-45 2 + {l-2B)xl 


2 5 - 8 B 

22C2B-1)x^ 


-3 2 F 1 


5-85 3 25-1 


2’ 2(1 - 25)’ 2’ 


-Xr 


+ Bxq 2 F 1 ( 


5-85 5 25 - 1 


2’ 2(1 -25)’ 2’ 


-Xr 


where 2-^1 (a, c; ^) is the Gauss hypergeometrie function defined by 

{a}n{b)n 


(6.4) 


2 Fi{a,b;c-,z) = ^ 


i =0 


(C)r 


n! 


with 


(a)o = 1, {a)n = a{a + 1)... (a + n - 1). 

The function 2 Fi{a, b; c; z) is analytic in the whole plane, except for a singularity at z = 1. 
Moreover, 2 Fi{a, b; c; 0) = 1 and 


25 i(a, b] c; z) = ^ ^ 2 Fi{a,b] a + b - c + 1]1 - z) 

1 (c — a)l (c — 0 ) 

+ (1 — zy “ ^ ^ ^ 2 .^i(c — a,c — b;c — a — b + 1;1 — z) 

T{a)T{b) 

for ^ G (—1,1). This yields that near z = 1 

p u. _ r(c)r(c -a-b) r(c)r(a + b-c) 

- r(c-a)r(c-6) + r(a)r(6) ' 

In the two hypergeometrie functions appearing in (6.4), the exponent of (1 — z) in ( 6 . 6 ) is 

- - J5^r<0, 5>|, 

2(1-25)1>0, 5<f.- 


(6.5) 


( 6 . 6 ) 


Hence, the first (resp. second) term in ( 6 . 6 ) is dominant when 5 < I (resp. 5 >!)■ We treat 
the two cases. 

The case 5 < |. For close to 1, the bracket part of (6.4) is close to 


,r(i)r(5iSI)) „ , r(|)r(5^) r(l)r(4^,) 






r(i)r(5S) ■ "“r(i)r(^ + i) r(i)r(ja=i) 

since r(a: + 1) = a:r(a:). We let x^ = — 5, with 5 > 0 small. Using that r(|) = |a/^, the 

integral part of Is in (6.3) is close to 

30r ^( 2 ( 1 - 2 !)) 2B-1 (n pp D _± 1 
2(5-1) (^ ~ ^0 7^ i 

30Fr(2()L-2l)) p^2 

4 r(^) 


5x2 +0(5-5o), 5o = l. 
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The coefficient is nonzero for 5 > 0 as soon as -Bq 7 ^ 1 (resp. Bq = 1) and (resp. 
— is not a negative integer, which is the case for B > This shows that Is grows as 

(2 + (1 — 2B)xl)'^'^^-^^'>, while 2 (/i + I 2 ) grows as (2 + (1 — 2B)xl)~^. Hence, [[I 3 ]] is 
dominant when B <l and 2 (Ji + I 2 ) + I 3 = 7 ^ 0 when 5 < |. 

The case -B > |. [[For close to 1, the bracket part of (6.4) has two parts J 3 and J 3 . 




4B-3 



This yields the corresponding parts J 3 and Jg for J3, considering that 6 = 2 + {1 — 2B)xl: 



Considering that 2 (i- 2 B) ^ 0)’ + -^ 2 ) + h = (1 + 0(5)) 7 ^ 0. ]] □ 
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